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Abstract: The MAGIC telescopes discovered very high energy (VHE, E>100 GeV) gamma-ray emission coming from 
the distant Flat Spectrum Radio Quasar (FSRQ) PKS 1222+21 (4C +21.35, z=0.432). It is the second most distant VHE 
gamma-ray source, with well measured redshift, detected until now. The observation was performed on 2010 June 17 
(MJD 55364.9) using the two 17 m diameter imaging Cherenkov telescopes on La Palma (Canary Islands, Spain). The 
MAGIC detection coincides with high energy MeV/GeV gamma-ray activity measured by the Large Area Telescope 
(LAT) on board the Fermi satellite. The averaged integral flux above 100 GeV is equivalent to 1 Crab Nebula flux. The 
VHE flux measured by MAGIC varies significantly within the 30 minutes of exposure implying a flux doubling time of 
about 10 minutes. The VHE and MeV/GeV spectra, corrected for the absorption by the extragalactic background light, 
can be described by a single power law with photon index 2.72 ± 0.34 between 3 GeV and 400 GeV, consistent with 
gamma-ray emission belonging to a single component in the jet. The absence of a spectral cutoff at 30-60 GeV (indeed, 
one finds a strict lower limit E c > 130 GeV) constrains the gamma-ray emission region to lie outside the broad line 
region, which would otherwise absorb the VHE gamma-rays. Together with the detected fast variability, this challenges 
present emission models from jets in FSRQs. 
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1 Introduction 

Blazars are active galactic nuclei hosting powerful rela- 
tivistic jets pointing toward the observer. They are char- 
acterized by strong non-thermal emission extending across 
the entire electromagnetic spectrum, from radio up to 7- 
rays. The blazar PKS 1222+21 (also knows as 4C 21.35) 
belongs the class of Flat Spectrum Radio Quasars (FSRQs) 
and is located at a redshift of z = 0.432 [1]. FSRQs dis- 
play luminous broad emission lines, often accompanied by 
a "blue bump" in the optical-UV region which is associ- 
ated with the direct emission from the accretion disk. This 
optical-UV emission is believed to be re-processed by the 
Broad Line Region (BLR) clouds, filling the BLR with a 
dense optical-UV radiation field which can interact with 
VHE 7-rays causing internal absorption in the VHE do- 
main [6]. PKS 1222+21 is the third FSRQ detected in VHE 
7-rays after PKS 1510-08 [2] and 3C 279 [3, 4] at the red- 
shifts of z = 0.36 and z — 0.536, respectively. It is thus 



the second most distant VHE source after 3C 279 with well 
measured redshift (the BL Lac 3C 66A with an estimated 
redshift z = 0.444 [5] would occupy this position of the 
ranking, but the measurement has a large uncertainty). 

In this proceeding, we will present the MAGIC discovery 
of this source which was detected during a 7-rays flare an- 
nounced by the FerrrufLAT collaboration and its physics 
implications. 

2 MAGIC discovery 

MAGIC consists of two 17 m diameter Imaging Atmo- 
spheric Cherenkov Telescopes (IACTs). The data were 
taken at zenith angles between 26° and 35°. Stereoscopic 
events, triggered by both telescopes, were analyzed in the 
MARS analysis framework [8]. Details on the analysis can 
be found in [23] whereas the performance of the MAGIC 
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Figure 1: Distribution of the squared angular distance (8 2 ) for 
events in the direction of 4C +21.35 (black points) and 
normalized off-source events (grey histogram). The signal is 
extracted in the 8 2 -region denoted by the vertical dashed line. 

stereo system will be discussed in detail in a forthcoming 
paper [10]. 

PKS 1222+21 was observed by MAGIC on June 17 (MJD 
55364) for 30 minutes as part of a Target of Opportunity 
program triggered by an increased flux in the Fermi en- 
ergy band [9]. During this detection by MAGIC the source 
was close to the brightest flare ever observed by the Fermi 
Large Area Telescope (LAT) [22]. The signal evaluation 
was performed using the 8 2 distribution (squared angular 
distance between the true and reconstructed source posi- 
tion), see Fig. 1. We got an excess of 190 7-like events 
(6 7/min.) above a background of 86 events, which corre- 
sponds to a statistical significance of 10.2 a using eq. 17 in 
[7]. The energy threshold of this analysis is s=s 70GeV. 

2.1 Very High Energy Spectrum 

The differential energy spectrum measured by MAGIC ex- 
tends from 70GeV to 400 GeV (Figure 2) and it is well- 
described by a simple power law of the form: 
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with a photon index r = 3.75 ±0.27 sla ,±0.2 syst and a nor- 
malization constant at 200 GeV of A^oo = (7.8 ± 1.2 sla ,± 
3.5 sysl ) x 10 _10 cm _2 s _1 TeV _1 , yielding an integral flux 
(4.6 ± 0.5) x lO-^cm^s- 1 (« 1 Crab Nebula flux) at 
E> 100 GeV. For energies higher than 400 GeV no signif- 
icant excess was found but two upper limits correspond- 
ing to 95% confidence level (C.L.) have been derived. The 
method used for the spectral reconstruction is the so-called 
"Tikhonov" unfolding algorithm [15], which takes into ac- 
count the finite energy resolution of the instrument and the 
biases in the energy reconstruction. The systematic un- 
certainty of the analysis (studied by using different cuts 
and different unfolding algorithms) is shown by the striped 
area. 
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Figure 2: Differential energy spectrum of PKS 1222+21 as mea- 
sured by MAGIC on 2010 June 17. Differential fluxes are shown 
as black points, upper limits (95% C.L.) as black arrows. The 
black line is the best fit to a power law. The black-striped area rep- 
resents the systematic uncertainties of the analysis. The absorp- 
tion corrected spectrum and upper limits using the EBL model of 
[11] are shown by the grey squares and arrows; the dashed grey 
line is the best fit power law. The grey shaded area illustrates the 
uncertainties due to differences in the EBL models of [12, 13, 14] 
and [3]. 



The deabsorbed spectrum is shown by the grey squares in 
Fig. 2, where the EBL model of [11] has been used. It is 
well fitted by a power law with an intrinsic photon index of 
r tott = 2.72 ±0.34 between 70 GeV and 400 GeV. Uncer- 
tainties caused by the differences between the EBL models 
are represented in Fig. 2 by the grey shaded area. The cor- 
responding spread is smaller than the systematic uncertain- 
ties of the MAGIC data analysis. 

Using the \ 2 difference method we studied the possibility 
of a cut-off in the energy spectrum. We tried to fit the spec- 
trum by a broken power law with different photon indexes 
and values for the cut-off. We concluded from this study 
that with the available statistics, at the 95% C.L. we cannot 
exclude the presence of a cut-off above 130 GeV for a pho- 
ton index 2.4 (the lowest possible value compatible with fit 
uncertainties and with the Ferm//LAT data, see Fig. 4) or 
above 180 GeV for a photon index 2.7. 

2.2 Variability in Very High Energy 7-rays 

Thanks to the strength of the signal even if the observa- 
tion time is as short as 30 minutes, a variability study is 
possible. In Fig. 3 the light curve binned in 6 minute long 
intervals is shown. It reveals a flux variation within the 30 
minutes of observation time. The hypothesis of a constant 
flux is rejected (x 2 /NDF = 28.3/4) with high confidence 
(probability < 1.1 x 10~ 5 ). The flux of the background 
events surviving the 7/hadron selection cuts is compatible 
with being constant and hence we can exclude a variation 
of the instrument performance during the observation. 
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Figure 3: PKS 1222+21 light curve above 100 GeV, in 6 minute 
bins (black filled circles). The observation was carried out on 
MJD 55364. The black solid line is a fit with an exponential func- 
tion and the black dotted line a fit with a linear function. The 
grey open squares denote the fluxes from the background events 
and the grey dashed line is a fit with a constant function to these 
points. 

To quantify the variability time scale we performed an ex- 
ponential fit (solid black line in Fig. 3). A linear fit is 
also acceptable but does not allow us to define a time scale 
unambiguously. For the exponential fit the doubling time 
of the flare is estimated as 8. 6+ 9 minutes. The derived 
timescale corresponds to the fastest time variation ever ob- 
served in an FSRQ in the VHE range and in any other en- 
ergy range [20], and is amongst the shortest measured from 
any TeV emitting source [21, 17]. 

3 Spectral Energy Distribution 

The high-energy Spectral Energy Distribution (SED) is 
shown in Figure 4 where the MeV/GeV energy range spec- 
trum from Fermi/LAT and the GeV/TeV spectrum mea- 
sured by MAGIC are combined. The source showed a 
significant flare in the Fermi band lasting ^3 days, with 
a peak flux on 2010 June 18 (MJD 55365) [22]. During 
the 30 minutes MAGIC observation there is a gap in the 
LAT exposure, so we analyzed the closest data available, a 
period of 2.5 hr (MJD 55364.867 to 55364.973) before and 
after the MAGIC observation. Given the short observation 
time (chosen in order to be as much contemporaneous as 
possible with the MAGIC data) there is no detection above 
2 GeV in the Fermi/LAT data, but an upper limit at the 95% 
C.L. in the energy range 2 — 6.3 GeV has been calculated 
together with the spectral points up to 2 GeV and combined 
with the MAGIC data in the SED shown in Fig. 4. 

If we extrapolate the intrinsic MAGIC spectrum to lower 
energies we can see that there is a potentially smooth 
connection between the Fermi/LAT and MAGIC extrapo- 
lated data and the photon index steepens from 1.9 in the 
Fermi/LAT range to 2.7 in the MAGIC range. These results 
agree with the analysis of larger temporal intervals during 




Figure 4: High energy SED of PKS 1222+21 during the flare 
of 2010 June 17 (MJD 55364.9), showing Fermi/LAT (squares) 
and MAGIC (circles) differential fluxes. A red bow tie in the 
MeV/GeV range represents the uncertainty of the likelihood fit 
to the Fermi/LAT data. The unfolded and deabsorbed spectral fit 
of the MAGIC data is also shown as a grey bow tie, extrapolated 
to lower and higher energies (dotted lines). A thick solid line 
(photon index T — 2.7) indicates a possible extrapolation of the 
MAGIC deabsorbed data to lower energies. The thick dashed line 
represents the EBL absorbed spectrum obtained from the extrap- 
olated intrinsic spectrum using the model of [11]. 

this flare and during the whole active period, in which the 
source spectrum is well described by a broken power law 
with an energy break between 1 and 3 GeV [22] . Further- 
more it is found that the high energy tail (E > 2 GeV) of the 
Fermi/LAT spectrum of 4C +21.35 extends up to 50 GeV, 
with a photon index in the range 2.4-2.8. 

4 Physical implications on the jet emission 
models 

As we have discussed, the almost simultaneous VHE and 
GeV spectra are consistent with a single power law with 
index ~ 2.7 ± 0.3 between 3 GeV and 400 GeV, without 
a strong intrinsic cutoff and with a smooth connection to 
the lower energy spectrum. This evidence suggests that the 
100 MeV - 400 GeV emission belongs to a unique compo- 
nent, peaking at « 2 — 3 GeV, produced in a single region 
of the jet. Considering the inverse Compton scattering on 
external photons and relativistic electrons in the jet as the 
emission process, as it is usually assumed, we have two 
possible scenarios. The emitting jet region could be in- 
side of the Broad Line Region (BLR), where the external 
photons field would be the UV photons from the BLR or 
we can assume this emitting region to be outside of the 
BLR where the photon field would be composed of the IR 
photons coming from the torus. There are two important 
effects that need to be taken into account: the decreased ef- 
ficiency of the IC scattering occuring in the Klein-Nishina 
(KN) regime and the absorption of 7-rays through pair pro- 
duction. 
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The energy above which the KN effects and 7-ray absorp- 
tion become important in the case of UV external photons 
from the BLR is ~ tens of GeV, and thus we would ex- 
pect a cut-off at these energies if the emitting jet region is 
inside the BLR. But if we consider as target photons the 
ones coming from the IR torus this effect would start to be 
appreciated only at much higher energies above ss 1 TeV. 
Since there is no evidences of a cut-off at low energies we 
can conclude that the emission should come from outside 
of the BLR, as has been proposed by the "far dissipation" 
scenarios [e.g. [19]]. 

Besides, the other important result of the MAGIC observa- 
tion is the fast variability, i var ^10 minutes which indi- 
cates an extremely compact emission region. This is diffi- 
cult to reconcile with the "far dissipation" scenarios if the 
emission takes place in the entire cross section of the jet 
since in this case the emitting region should be close to the 
central black hole and thus inside of the BLR. 

Some explanations have been already proposed in order to 
solve this kind of incongruities, invoking the presence of 
very compact emission regions embedded within the large 
scale jet [18, 24] or the possibility of a very strong jet rec- 
ollimation [e.g. [16]]. 

In conclusion the MAGIC observation of VHE emission 
from the FSRQ PKS 1222+21 sets severe constraints on 
the emission models of blazar jets. 
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